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A Sampling Method for Determining the Lowest 
Eigenvalue and the Principal Eigenfunction of 
Schrodinger s Equation 1 

By M. D. Donsker and M. Kac 

This is a preliminary report on a sampling method for finding the lowest eigenvalue 
and the corresponding eigenfunction of Schrodinger's equation. The theory underlying 
the procedure is explained and two cases treated numerically. Although the initial results 
are encouraging, more experimentation will be needed to fully test the practicality of the 
method. 



I. Introduction 

Certain problems leading to complicated partial 
or integro-diff erential equations have recently been 
approached and some actually solved by utilizing 
various probability techniques and sampling 
methods. Collectively these methods have be- 
come known as the " Monte Carlo" method. 

The problems to which Monte Carlo techniques 
have been applied seem to be divided into two 
types. Typical of the first type is the problem of 
neutrons diffusing in material media in which the 
particles are subjected not only to certain deter- 
ministic influences but to random influences as 
well. In such a problem, the Monte Carlo ap- 
proach consists in permitting a " particle" to play 
a game of chance, the rules of the game being such 
that the actual deterministic and random features 
of the physical process are step by step exactly 
imitated by the game. By considering very large 
numbers of particles, one can answer such ques- 
tions as the distribution of the particles at the end 
of a certain period of time, the number of particles 
to escape through a shield of specified thickness, 
etc. One important characteristic of the preced- 
ing approach is that the functional equation de- 
scribing the diffusion process is bypassed com- 
pletely, the probability model used being derived 
from the process itself. 

A more sophisticated application of Monte 

1 The preparation of this paper was sponsored (in part) by the Office of 
Naval Research. 



Carlo methods is to the problem of finding a 
probability model or game whose solution is re- 
lated to the solution of a partial differential 
equation, or, as in the present paper, to determine 
the least eigenvalue of a differential operator by 
means of a sampling process. As an example of 
how the latter problem might be attacked, we 
quote from a paper of Metropolis and Ulam: 2 

"For example, as suggested by Fermi, the time 
independent Schrodinger equation 

A<l>(v,y,z) = (\—V)<l>(z,y,z), 



could be studied as follows, 
dependence by considering 



Reintroduce time 



u(x,y 1 z i t)=<j>{x,y,z)e-*) 
then, u will obey the equation 



dt 



=Au—Vu. 



This last equation can be interpreted, however, as 
describing the behavior of a system of particles 
each of which performs a random walk, i. e., 
diffuses isotropically and at the same time is 
subject to multiplication, which is determined by 
the value of the point function V. If the solution 
of the latter equation corresponds to a spatial 

2 N. Metropolis and S. Ulam, The Monte Carlo method, J. Am. Stat. 
Assn. 44, 247, pp. 335 to 341 (Sept. 1949). 
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mode multiplying exponentially in time, the ex- 
amination of the spatial part will give the desired 
<f>(x,y,z) — corresponding to the lowest 'eigenvalue' 

The main purpose of the present paper is to 
present an alternative method for finding the low- 
est eigenvalue and corresponding eigenfunction of 
Schrodinger's equation. The chief difference be- 
tween the two approaches is that ours involves 
only a random walk eliminating entirely the 
multiplicative process. This alteration in the 
model seems to simplify the numerical aspects of 
the problem, especially if punched card equipment 
is to be used. Apart from the possible numerical 
simplification, the method is based on a mathemat- 
ical theory that in itself is of some interest. 

II. Mathematical Theory- 
Let X u X 2 , X s , ... be independent identically 
distributed random variables each having mean 
and standard deviation 1 and let S k =Xi + X 2 + 
. . . -\-X k . Under certain general assumptions 
on V(x), the most severe of which is that V(x) be 
nonnegative, it can be shown 4 that the limiting 
distribution function <r(a,t) of the random variable 



n k<nt \^n/ 



is such that 



"e- a - st d a <r(a,t)dt= J x//(x)dx, 



(1) 



(2) 



where \p(x) is the fundamental solution of the 
differential equation 



1 dV 



%r-( 8 + V (x))4,=0, 



(3) 



2dx 

subject to the conditions 

^(x)—>0 x-^±<*> 

\+'(x)\<M x^O 

*'-(+0)-*'(-0) = --2. 

The fundamental solution \p(x) of (3) is expres- 
sible in terms of the normalized eigenfunction 

3 We have since learned that Dr. G. W. King has made similar calcula- 
tions using a closely related method. 

< M. Kac, On distributions of certain Wiener functionals, Trans. Am. 
Math. Soc. 65, 1 to 13 (1949). 



{(^0*0} a nd eigenvalues \j of the one-dimensional 
Schrodinger eigenvalue problem 5 



1 y/2 / 

2^-V(^(x)= -M X>0 (4) 

W=S^. (5) 



as 



Thus from (5 and 2) 

J CO /»co /»co 

_ f{x)dx= " e-"- s 'd a a{a,t)dt 

MO) f " *&)<** 

Inverting (6) with respect to s we get, 

/» 00 f* CO 

«-«rf a cr(a,0=S«- M ^(O) M*)dx, (7) 

JO J-oo 

and therefore we obtain the expression for X x , 

1 f °° 

Xx^lim— T log e- a d a <r(a,t). (8) 

*->«, t Jo 

If in (7) we neglect all terms in the expansion 
but the first, we obtain 

log I e~ a d a <T(a,t) ~logWi{0)\ rl/iffldxt—Xit, 
or 



Xi 



1 f ro 

-log e- a d a a(a,t). 

(9) 



Thus, if, by choosing a finite t, we attempt to cal- 
culate Xi from (8), we have two sources of error. 
The first, usually a small source of error, is from 
the exponentials neglected in the expansion (7). 
The second, and more important, is from neglect- 
ing the term (l/i)log\ypi(0)\_yi(x)dxl- This 

latter source of error is especially significant, since, 
as will be apparent shortly, it is impractical from 
other points of view to take t very large. All of 
this difficulty may be obviated by considering (7) 
for two distinct values of t, say t x and t 2 , then, if 



« From here on all the steps are formal. In all cases of physical interest 
they can be justified rigorously. 
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we neglect the exponentials after the first as before, 
we get on dividing 

K4V do) 



Xi' 



t 2 -t x 



e a d a <j{a, 
Jo 



U) 



The Monte Carlo process consists in the calcula- 
tion of a(a,ti) and <r(a,t 2 ) by a sampling process. 

If instead of a(a,t) we consider the limiting dis- 
tribution az(a,t) of the random variable 

1 2 vU+^f), (ii) 

n k<nt \ -yjn/ 

then <rt(a,t) also satisfies (2 and 3), but now the 
condition \p' ( + 0) — \p' ( — 0) = —2 is replaced by 
^'(£+) — ^' (£—) = - 2. Therefore, repeating steps 

(4 to 7), we get 

f"€-«rf^(a,0«S«-^i({) r°° **(*)<&• (12) 

J o j J - °° 



Thus 



*i(0) 



lim -^ 

Jo 



(13) 



so that we are also able to calculate the principal 
eigenf unction. 

The extension of the preceding method to multi- 
dimensional Schrodinger equations is immediate. 
It is in these cases that the method will probably 
prove to be most useful since, unlike the standard 
variational procedures, the extension to several 
dimensions seems to cause comparatively little 
difficulty. For illustrative purposes we will con- 
sider Schrodinger's. equation in three dimensions. 
Here we must consider three independent 
sequences 

Xi, X 2 , X$, . . . , 

Yl, JL2y Yzy . . . , 

Zi> Z 2 , Z 3 , . . . , 

of independent identically distributed random 
variables each having mean and standard devia- 
tion 1. Let S xk =X l +X 2 + - • • +X k) and S ykf 
S zk have the obvious meanings. We consider 
the limiting distributions <r(a,t) and a(a, {, rj, f, t) of 
the random variables 



and 



n k<nt V-y'Ti y/n -yln/ 
w k<nt \ \n ^n \ a ' 



respectively. 

In exactly the same way as in the one dimen- 
sional case, we arrive at 

f* oo 

l I er<*d a <j{a,t { ) 

Xl^T— 7-lOg "~oo ~, 

e-°d a <r(a,t 2 ) 



and 



*, (0,0,0) ". 



j o e-«d a <r{a,$,r,,S,t) 

T 

Jo 



e a d a a(a,t) 



So far, the theory was carried out under the 
assumption that the potential function V was 
nonnegative. In most cases of physical interest 
this is not so. For the hydrogen atom, for 
instance 

const. 



V(x,y,z)- 



V* 2 -f-2/ 2 +z 2 



However, the modification is easy, although we 
have not as yet cleared up all the points of math- 
ematical rigor. 

The formula for the lowest eigenvalue now 
becomes 

l e-<*d a (j{a,U) 

\_J- a daa(a,t 2 ) 

and a corresponding modification needs to be 
made in the formula for the principal eigenf mic- 
tion. We have not yet tested numerically any 
case with a negative potential function, but we 
hope to be able to report on this in the near future. 

III. Numerical Examples and Discussion 

The Monte Carlo procedure used here consists 
in the calculation of the distribution function 
cr(a,t) by a sampling process; the principal eigen- 
value is then calculated from (10). For the pur- 
poses of numerical illustration we consider two 
examples, V(x)=x 2 and T 7 (#) = |a;|. In both of 
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these cases the eigenvalues are known, and hence 
we have a check on the accuracy of our procedure. 
In order to calculate <r(a,t), say when V(x)=x 2 , 
we see from (1) that we must consider the limiting 
distribution as n-* °° of 



^=5, £2=3.75, and 100 samples were used, 
table 1 is tabulated 



In 



aS& ! 



n*k<Zt 



(14) 



This means from our point of view that we 
must (a) choose a distribution for the X's, (b) 
choose a sufficiently large n, (c) select an appro- 
priate t, (4) calculate for nt Js the normalized 
sum (14), (5) repeat (4) many times, so that the 
empirical distribution may be obtained from these 
many samples. 

Although, under the conditions mentioned 
previously, the distribution function a(a,t) is 
independent of the distribution of the X's, the 
actual numerical calculation of <r(a,t) is expedited 
by choosing the distribution of the X's to be the 
Bernoulli distribution, i. e., 



P{X=l)=P(X=-l)=2 



The sequence of random variables Xi, X 2 , X z , 
... is then a sequence of + l's and — l's, such 
as might be obtained in coin tossing. This is 
conveniently and rapidly achieved on a calculat- 
ing machine by considering sequences of random 
digits, counting even digits +1 and odd digits — 1. 
The value of n to be used must be large enough 
so that the empirical distribution function cal- 
culated is close to the theoretical limiting dis- 
tribution function <r(a,t). From (10) we see that 
the two values of t, t h and t 2 , to be used in the 
calculation of X must be large enough so that the 
exponential terms neglected are sufficiently small. 
However, since the sample size is nt, the desire 
to make both n and t large must be tempered by 
practical considerations. The number of samples 
to be used must be large enough so that the 
empirical distribution adequately represents a(a,t). 
Before discussing these points in more detail we 
consider an actual numerical computation. The 
following data for V(x)=x 2 and F(a?) = |aj| were 
calculated from a certain set of random digits 6 
on the IBM Electronic Calculating Punch, Type 
604. For both x 2 and \x\ n was selected to be 400, 



Av 



1 
1500 



I 2000 



1500 

B, = 160000 S S **' W= 



2000 
)0/c = 



SW 



Since each rand random number card contains 
50 random digits, and 2,000 digits are needed to 
form one sample, a set of 4,000 rand cards was 
sufficient for this experiment. It takes 20 min- 
utes on the IBM 604 to calculate A x and A 2 for 
one sample (similarly for B x and B 2 ), so that it 
takes approximately 35 hours to secure the follow- 
ing data. 

Table 1. Experimental data for lowest eigenvalue 



• This 'set of random digits was prepared by the RAND Corporation, 
Santa Monica, Calii. 





Ax 


A* 


Bx 


B 2 


1 


1.332 


2.101 


0.738 


1.290 


2 


3.808 


5.707 


5.295 


8.502 


3 


2.795 


5.545 


3.460 


9.610 


4 


8.723 


13. 917 


30. 722 


52. 616 


5 


4.169 


4.766 


6.123 


6.602 


6 


1.195 


2.136 


.598 


1.453 


7 


6.674 


12.112 


15. 258 


39.058 


8 


4.103 


5.242 


5.826 


7.004 


9 


5.751 


8.840 


10. 440 


18. 143 


10 


4.250 


4.981 


6.069 


6.679 


11 


2.909 


5.643 


3.844 


10.044 


12 


2.834 


3.416 


3.602 


4.097 


13 


1.888 


2.194 


1.595 


1.700 


14 


2.022 


2.337 


1.510 


1.638 


15 


1.680 


3.908 


1.184 


5.289 


16 


7.700 


12.712 


24. 769 


44. 936 


17 


3.228 


4.973 


3.563 


6.145 


18 


1.844 


2.654 


1.523 


2.183 


19 


2.376 


5.275 


2.017 


8.854 


20 


4.533 


5.640 


6.380 


7.616 


21 


4.209 


7.543 


7.735 


16. 703 


22 


3.847 


5.590 


6.666 


9.204 


23 


5.206 


6.314 


9.148 


10. 600 


24 


1.962 


2.698 


1.859 


2.923 


25 


3.523 


4.807 


5.084 


6. 482 


26 


3.605 


5.397 


5.014 


7.680 


27 


9.168 


13. 692 


27. 973 


44. 461 


28 


5.625 


7.898 


11. 528 


16. 010 


29 


3.215 


4.152 


3.362 


4.170 


30 


3.293 


3.972 


4.136 


4.677 


31 


2.710 


3.683 


3.006 


4.025 


32 


3.638 


7.386 


6.114 


17. 413 


33 


7.928 


11. 567 


25.084 


36.001 


34 


3.475 


3.846 


4.396 


4.550 


35 


1.958 


2.594 


1.457 


1.948 


36 


5.356 


7.216 


[ 8. 888 


11. 781 


37 


3.297 


6.395 


3.800 


11.555 


38 


1.883 


2.682 


1.387 


1.995 


39 


4.849 


10.020 


8.396 


30. 233 


40 


2.282 


3.067 


2.258 


2.882 
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Table 1. Experimental data for lowest eigenvalue — Con. 





A\ 


At 


Bi 


Bi 


41 


4.253 


6.957 


6.744 


13.502 


42 


7.559 


8.834 


17. 143 


18. 702 


43 


1.122 


2.587 


.500 


2.634 


44 


1.433 


3.602 


.869 


4.855 


45 


2.984 


3.477 


2.851 


3.161 


46 


1.677 


3.158 


1.025 


3.341 


47 


5.145 


9.697 


11. 941 


28 826 


48 


4.982 


9.140 


9.201 


23. 092 


49 


1.675 


2.152 


1.129 


1.472 


50 


1.580 


3.375 


1.147 


3.807 


51 


2.606 


3.621 


2.254 


3.247 


52 


4.303 


5.883 


7.987 


10. 104 


53 


7.553 


12. 092 


19. 632 


36. 427 


54 


3.138 


5.703 


3.968 


9.325 


55 


3.343 


4.327 


4.264 


5.380 


56 


1.423 


2.315 


0.770 


1.533 


57 


7.996 


11. 458 


21.405 


31. 585 


58 


1.916 


2.744 


1.253 


1.959 


59 


1.490 


2.161 


0.789 


1.404 


60 


3.536 


5.382 


4.125 


7.461 


61 


6.278 


11.060 


16. 187 


34. 779 


62 


1.174 


1.709 


0.573 


0.919 


63 


5.272 


9.919 


10.463 


27. 825 


64 


3.576 


5.277 


6.021 


8.707 


65 


2.686 


5.081 


3.577 


8.307 


66 


6.166 


8.785 


14.056 


19. 627 


67 


4.412 


5.548 


6.095 


7.182 


68 


1.732 


4.664 


1.170 


8.720 


69 


2.730 


4.970 


2.440 


6.639 


70 


1.331 


2.772 


0.765 


2.642 


71 


3.668 


5.318 


5. 730 


8.001 


72 


2.948 


3.801 


3.248 


4.104 


73 


5.279 


8.931 


12. 455 


23. 513 


74 


11.284 


15. 784 


41.684 


58. 525 


75 


2.183 


2.787 


2.086 


2.455 


76 


4.298 


7.427 


9. 139 


17. 157 


77 


1.412 


2.893 


0.894 


3.490 


78 


2.002 


2.567 


1.602 


2.024 


79 


2.619 


3.780 


2.981 


4.205 


80 


2.695 


4.616 


2.974 


6.077 


81 


4. 706 


6. 672 


8.330 


11.487 


82 


5.517 


7.981 


9.846 


14. 805 


83 


1.755 


2.099 


1.066 


1.207 


84 


2.084 


2.876 


1.735 


2.504 


85 


1.574 


2.743 


0.993 


2.256 


86 


5.818 


10. 444 


13. 322 


30. 522 


87 


4.387 


7.844 


9.037 


18. 695 


88 


3.666 


5.244 


4.416 


6.541 


89 


7.111 


13.311 


19. 834 


50. 822 


90 


6.280 


8.548 


13. 403 


17. 587 


91 


4.872 


5.836 


9.230 


10. 023 


92 


3.035 


3.323 


2.865 


2.973 


93 


3.608 


4.697 


4.256 


5.357 


94 


3.977 


4.964 


5.215 


6.242 


95 


3.614 


7.867 


5.757 


20.434 


96 


9.193 


11. 183 


27. 789 


31.075 


97 


4.108 


5.282 


5.918 


7.374 


98 


4.355 


5.744 


7.774 


10. 223 


99 


1.898 


2.661 


1.507 


2.013 


100 


9.683 


17. 154 


35. 803 


80. 543 



In order to calculate Xi from (10) we need the values 

of e'^ao-ia^) and e~ a d a <j(a,Q. Both of these 

integrals were calculated numerically from the 
data in table 1 by adding the exponentials of the 
entries in appropriate columns. For the case 
V(x) = \x\, the values of the integrals were obtained 

100 100 

from 1/100 2>"" A " and 1/100 S f*" and for V{x) 

-•=1 ;=1 

100 100 

=x 2 from 1/100 X! e~ Bli and 1/100 2 e ~ B2i - In tIie 
j=i j=i 

case V(x) = |x| the true lowest eigenvalue to two 
places is 0.81, and in the case V(x) =x 2 it is-y/2/2= 
0.71. 





First 50 
samples 


Second 50 
samples 


All 100 
samples 


1*1 

X 2 


0. 83 
.80 


0. 79 
. 69 


0. 81 
.75 



An interesting feature of the data is that column 
B 2 , for example, may be thought of as arising from 
^ = 2000, t=l) n = 1000, *=2, etc., as well as from 
n=400, 2=5. The larger we take n, the closer 
will be the empirical distribution to <r(a,f). The 

value of e~ a d a (r(a,t) in the case t=l, and V(x) 

=x 2 can be calculated exactly and, to three places, 
is 0.678. Calculating this integral from column 
B 2 with n = 2000, t=l, we get 0.685. 

Instead of using Bernoulli distributed variables, 
one might use other distributions. One definite 
advantage of Bernoulli distributed variables is 
that the computation utilizes only the crudest 
properties of the random digits, i. e., whether they 
are even or odd. One possible advantage for cer- 
tain other distributions is that n might not have 
to be taken so large. In particular this should 
be true if we use Gaussian distributed variables. 
rand Gaussian deviates were used in constructing 
table 2. Here we chose /i = 3.75, fe=5 as before, 
but now n=10Q. This means samples of size 500 
instead of 2000 and therefore a total machine 
computation time of 9 hours for 100 samples. 

Envisaging the possibility of calculating the 
second eigenvalue, we considered the quantities 

375 500 

10- 4 S 0&+5) 2 and 10-" ZJ (S*+5) 2 , (15) 

k=l lt=l 
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Table 2. Experimental data for second eigenvalue 



Table 2. Experimental data for second eigenvalue — Con. 



1 


c, 


c 2 


61 


C 


Ci 


3.963 


13. 997 


49. 929 


84. 307 


2 


2. 348* 


2.721 


62 


3.095 


6.220 


3 


1. 889* 


2.231 


63 


2. 585* 


3.867* 


4 


3.817* 


4. 372* 


64 


1.318* 


1. 906* 


5 


2.395 


2. 802* 


65 


2. 241* 


9. 102* 


6 


12. 467 


27. 683 


66 


9.737 


12. 185 


7 


1. 068* 


1. 115* 


67 


9. 738* 


11.347* 


8 


0. 578* 


0. 967* 


68 


4.604 


5.589 


9 


1. 504* 


1.760 


69 


20. 364 


29. 707 


10 


4.249 


9.787 


70 


1.097* 


1. 895* 


11 


1. 751* 


3. 316* 


71 


5.891 


13.391 


12 


9.922 


13. 975 


72 


9. 019* 


12.815* 


13 


5.680 


6.649 


73 


1.843 


6.280 


14 


3.348 


4.684 


74 


7.961 


11.431 


15 


13. 431 


24. 769 


75 


1.109 


1.926 


16 


1.473* 


2.615* 


76 


2. 417* 


3. 037 


17 


45. 262 


74. 304 


77 


2. 402* 


2. 797* 


18 


1.157* 


3. 109 


78 


2.043* 


4.214* 


19 


0. 906* 


3. 644* 


79 


0. 582* 


0. 752 


20 


4.601 


12. 641 


80 


2.429 


4.381 


21 


8. 356* 


13. 084* 


81 


8.584 


13. 054 


22 


0. 308* 


1. 874* 


82 


2. 793* 


4. 659* 


23 


2. Ill* 


4. 428* 


83 


1. 557 


5.023 


24 


2. 035* 


2.726 


84 


1. 006* 


1. 677* 


25 


1. 625* 


3.898 


85 


5.327* 


15.531* 


26 


1. 879* 


5. 713 


86 


2. 058 


7.307 


27 


1.333 


8.961 


87 


12. 481 


12.803* 


28 


1.769 


2.319* 


88 


0. 925* 


2.403 


29 


1.906* 


2. 417* 


89 


15.519 


19. 879 


30 


7.522 


7. 979* 


90 


1. 456* 


4. 009* 


31 


14. 934 


15. 493* 


91 


5.244 


10. 047 


32 


4. 282* 


13. 747* 


92 


2. 497 


2.931 


33 


1.356* 


2.718 


93 


16. 700 


17.382 


34 


4.448 


8.329 


94 


4.471* 


8. 156* 


35 


1. 306* 


2. 112* 


95 


2. 139* 


2. 580* 


36 


5. 465* 


12. 396* 


96 


1. 196* 


2. 273* 


37 


2.918 


9.006 


97 


23. 014 


28. 420 


38 


35. 821 


56. 626 


98 


6. 959* 


7. 022* 


39 


9.814* 


11. 707* 


99 


3.374* 


8. 651* 


40 


1.811 


4.628 


100 


1. 600* 


2. 250* 


41 


3.374 


6.460 


101 


4. 626* 


7. 794* 


42 


1.633 


6.578 


102 


5. 195* 


10.215* 


43 


6.976 


8.705 


103 


6. 837* 


12. 160* 


44 


7. 595 


13. 096 


104 


1. 913 


4.412* 


45 


2.336 


2.862 


105 


2. 276* 


2.770 


46 


1. 931* 


2. 579* 


106 


1. 935* 


3. 983* 


47 


4. 507* 


8. 294* 


107 


3.560 


5.808 


48 


8.581* 


11.301* 


108 


2. 857* 


7.384* 


49 


4. 698* 


8. 220* 


109 


3.399* 


4. 484* 


50 


10. 154 


16. 884 














110 


6.917* 


27. 107* 


51 


8.828 


15. 387 


111 


38. 661 


56. 479 


52 


0. 997* 


1. 604* 


112 


1.891 


2. 000* 


53 


8.582 


12. 838 


113 


8.506 


9. 646* 


54 


2. 730* 


3. 249* 


114 


4.209 


10. 461 


55 


10. 471 


13. 140 


115 


33. 194 


44. 247 


56 


1. 456* 


1. 845* 


116 


0. 869* 


1.289 


57 


8. 230* 


11. 924* 


117 


3.946 


15. 204 


58 


5. 515* 


8. 146* 


118 


1. 680* 


2. 995* 


59 


2. 202* 


3. 182* 


119 


8. 178* 


9.056 


60 


2. 129* 


9. 650* 


120 


27. 621 


51. 650 



121 


C 


c 2 


161 


c, 


C 2 


8. 444* 


15. 836* 


3.546 


4.040* 


122 


13. 254 


16. 467 


162 


1. 238* 


2. 304* 


123 


1. 755* 


5. 246* 


163 


0. 975* 


1.283 


124 


36. 902 


70. 134 


164 


1.980 


3. 659* 


125 


7. 567* 


18. 553* 


165 


14. 620* 


17. 565* 


126 


.14. 725 


38. 551 


166 


2. 378* 


2.865 


127 


0. 983* 


2.572 


167 


2.279 


6.279 


128 


24. 529 


31. 660 


168 


10. 846 


13. 754 


129 


2. 042* 


2. 385* 


169 


1.341* 


3. 673* 


130 


3. 587* 


9. 858* 


170 


1.861* 


2. 729* 


131 


1.335* 


1.755 


171 


28. 588* 


53. 358* 


132 


2. 142* 


3. 801* 


172 


2.943 


3.087* 


133 


16.361 


27. 759 


173 


1.694 


2.568 


134 


8.193 


9.301 


174 


1. 729 


2.516* 


135 


2. 719* 


5. 141* 


175 


0. 844* 


1. 449* 


136 


4. 767* 


12. 737* 


176 


3.441* 


5. 380* 


137 


42. 230 


108. 038 


177 


0.712 


0. 838* 


138 


41. 236 


78. 423 


178 


1. 222* 


3.469* 


139 


8.243 


9. 911* 


179 


20. 852* 


42.263* 


140 


31. 522 


36. 393 


180 


7. 575* 


16. 552* 


141 


2.600 


2. 881* 


181 


2. 042* 


3.058 


142 


1.932 


2. 055* 


182 


6.797 


7. 343* 


143 


8.562 


10. 813 


183 


1.536 


4.499 


144 


1.421* 


3. 652* 


184 


1. 235 


1.472* 


145 


2. 669* 


3. 889* 


185 


3. 541* 


4.272 


146 


1.125* 


2.071 


186 


5. 394* 


6. 149* 


147 


3. 439* 


5. 824* 


187 


18. 589 


28. 787 


148 


1. 269* 


5. 651 


188 


8.372 


12. 628 


149 


7. 949 


14. 252 


189 


13. 470* 


23. 775* 


150 


2.808* 


4. 434* 


190 


1. 196* 


2.851* 


151 


6.413 


12. 573 


191 


3.507 


7.712 


152 


5. 936* 


10. 432* 


192 


21. 489 


43. 372 


153 


9.788 


18. 093 


193 


1. 038* 


1. 109* 


154 


1.251 


2.707 


194 


3. 591* 


3.864 


155 


34. 799 


63. 494 


195 


2. 295* 


2. 529* 


156 


0.841* 


1. 742* 


196 


1. 044* 


2.320 


157 


1.461* 


1. 615* 


197 


4. 865* 


6. 759* 


158 


7.220 


21. 125 


198 


1.184 


11.852 


159 


1. 275* 


2. 359* 


199 


1. 228* 


1. 969* 


160 


9.355 


19. 565 


200 


9. 346* 


22. 204* 



which correspond to taking £=0.5 in (11). This 
should not make any difference in the calculation 
of the lowest eigenvalue, and consequently table 2 
can be utilized in the same way as table 1. It 
should however be borne in mind that columns 
C\ and C 2 of table 2 represent experimental values 
of the quantities (18) with the S k 's being sums of 
Gaussian deviates. Asterisks on the entries of 
table 2 indicate that £375 or S 500 was negative. 
Although this information is unnecessary for the 
purpose of calculations of the lowest eigenvalue, 
it is used in the calculation of the next eigenvalue. 
How this can be done is explained briefly in 
section IV. 
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Using the data in table 2 and again (10) we 
computed the following approximations to Ai 
(actual value = 0.71) 



Xi 


First 50 
samples 


Second 50 
samples 


Third 50 
samples 


Fourth 50 
samples 


All 200 
samples 


0.82 


0. 72 


0. 83 


0. 64 


0. 74 



For both V(x) = \x\ and V(x)=x 2 , all the 
eigenvalues are known so that in the two illus- 
trative examples above it was easy to choose 
appropriate values for t x and t 2 . The proper 
value for n and the appropriate number of samples 
was determined experimentally. In general when 
all eigenvalues are unknown the following "rule 
of thumb" procedure seems to be feasible. Hav- 
ing first made a guess at an appropriate n and t, 
and a certain number of samples, compute 1 Xi. 
Repeat the computation now for the same n, 
the same number of samples and successively 
larger fs until the calculated values of Xi become 
stable to the desired number of places. If they 
do not become stable, the number of samples 
must be increased. Keeping two values of t, for 
which the calculated X's had the stable value, 
increase n and see if the calculated value of X 
changes. If not, n is sufficiently large. If it 
does change, increase n until a new stability 
appears. This stable value is then the appro- 
priate approximation to Xi. The mere fact that 
stability is observed means the number of samples 
is sufficiently large. 

The Second Eigenvalue 

If the principal eigenfunction is even, then it is 
possible to extend the theory of section II in 
such a way that the calculation of the second 
eigenvalue becomes feasible. Without going into 
any details we just state the pertinent result: 
Let 



Et(& = l\mE)e n k " nt V Vn7 |&,>o}' 

71— >oo \ 



and 



Er(Q = \\mE)e n k ^ nt V Vn/ |&,<o}' (17) 



the mathematical expectations on the right hand 
sides being conditional expectations under the 
conditions S nt ^>0 and aS^<C0, respectively. 
We then have 

\ { £"+(!) -Ef (© J =± e-^Mt) f\(x)fj{x)dx, 

(18) 
where 

( + l,z>0 

P(x) = ] 

(-l,z<0. 

If t l and t 2 are sufficiently large, we get 



1 



ff+(S)-ff,7(S) 

h-tr B EZ{i)-ErM) 



log 



(19) 



From the discussion of section III it should be clear 
how one applies (19) when the data of table 2 are 
available. One must see to it that £ is so chosen 
that ^ 2 (£)^0, otherwise a higher eigenvalue may 



have been calculated. 
we obtain 



From the data of table 2 



1.1, 



whereas the exact value is 3/2 V2 = 2.12. The poor 
agreement could have been expected in view of low 
accuracy in the calculation of \. 

In conclusion, we thank E. C. Yowell, of the 
National Bureau of Standards, for wiring the 
boards and for excellent supervision of all the 
punched card work. 

Los Angeles, September 12, 1949. 
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